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First order phase transitions in finite systems can be defined through the bimodal- 
ity of the distribution of the order parameter. This definition is equivalent to the 
one based on the inverted curvature of the thermodynamic potential. Moreover we 
show that it is in a one to one correspondence with the Yang Lee theorem in the 
thermodynamic limit. Bimodality is a necessary and sufficient condition for zeroes of 
the partition sum in the control intensive variable complex plane to be distributed on 
a line perpendicular to the real axis with a uniform density, scaling like the number 
of particles. 
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Phase transitions are universal phenomena which have been theoretically understood 
at the thermodynamic limit of infinite systems as anomalies in the associated equation of 
state (EoS). They have been classified according to the degree of non-analyticity of the 
thermodynamic potential at the transition point. For finite systems, it is generally believed 
that phase transitions cannot be denned since thermodynamic potentials of a finite system 
in a finite volume are analytic functions. This situation can be thought as unsatisfactory 
since on one hand the thermodynamic limit does not exist in nature, and on the other 
hand many experimental efforts are devoted to understand the link between mesoscopic and 
macroscopic systems, especially for spectacular properties like phase transitions |T|. 

A classification scheme valid for finite systems has been proposed by Grossmann |2j 
using the distribution of zeroes of the canonical partition sum in the complex temperature 
plane. Alternatively, it has been claimed that first order phase transitions in finite systems 
can be related to a negative microcanonical heat capacity |]3], |4[] or more generally to an 
inverted curvature of the thermodynamic potential as a function of an observable which 
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can then be seen as an order parameter ||. These anomalies can also be connected to 
the general topology of the potential energy surface ||. These very different approaches 
share the common viewpoint that at the mesoscopic or microscopic scale phase transitions 
are not smeared to become loosely defined smooth state changes, but can be analyzed in 
a rigorous way without explicitly addressing the thermodynamic limit. If this is certainly 
important to progress in the experimental study of phase transitions in finite systems, it is 
also true that the connection between these mesoscopic phenomena and the thermodynamics 
of the macroscopic world is not clear. This is especially annoying since the definitions of 
phase transitions in finite systems are not in general equivalent and different conclusions 
can be drawn out of the observation of a physical phenomenon depending on the theoretical 
framework used. In particular it can be shown |7j that any generic channel opening may lead 
to a convexity anomaly of the microcanonical entropy which does not necessarily survive to 
the bulk. In the same way the distribution of zeroes close to the real axis may suggest a first 
order phase transition even for systems where the transition temperature diverges in the 
thermodynamic limit ||. In order to have a coherent picture of the physical meaning of a 
phase transition in a finite system it is then important to analyze in detail the conditions of 
persistency of the signals towards the bulk and to make a clear bridge between the various 
definitions. 

We have recently proposed the possible bimodality of the probability distribution of 
observable quantities as a connection between these ideas ||. In this letter we further 
demonstrate that this definition is consistent with the Yang Lee theorem i.e. with the 
standard definition of first order phase transitions in the thermodynamic limit. Since our 
definition is a generalization of the previous ones @, |3], |J, the conditions of compatibility of 
the different signals and criteria proposed to identify first order phase transitions in finite 
systems can be explicitly worked out. 

Let us consider that for each event i we can measure a set of K independent observ- 
ables, b k , which form a space containing one possible order parameter. We can sort events 
according to the results of the measurements b^={b^} and thus define a probability dis- 
tribution P(b). In the absence of a phase transition logP(fe) is expected to be concave. An 
abnormal (e.g. bimodal) behavior of P{b) or a convexity anomaly signals a first order phase 
transition M. 

This definition does not imply an a priori knowledge of the thermodynamic potential, 
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i.e. of the properties of the specific equilibrium reached in the experimental situation, but 
only requires the measurement of the event distribution. Moreover, it does not assume the 
a-priori knowledge of an order parameter; the order parameter can be extracted from the 
data themselves by looking at the best direction in the observable space to split the event 
cloud into two components. This topological definition is therefore easy to experimentally 
implement, and in fact it has already allowed to successfully recognize phase transitions in 
metallic clusters and in the multifragmentation of atomic nuclei fl2| . 



Moreover this definition has a wider application domain than the previous ones 0, [3], |J ; 
it can be used out of canonical or microcanonical equilibria, for open systems || and non- 



extensive statistics such as Tsallis equilibria JT3[ . In particular it can be viewed as an exten- 



sion of the definition based on the convexity anomaly of the microcanonical entropy H |J or 
more generally of any thermodynamical Gibbs potential |J . Let us first discuss Boltzmann- 
Gibbs equilibrium which is obtained by maximizing the Shannon information entropy under 
the constraints of the various observables b k known in average. Then P(b) can be written 

as 

K 

log P x (b) = log\V(b) -J2 X kb k - \ogZ x (1) 

k=l 

where Z x is the partition sum of this intensive ensemble controlled by the Lagrange multipli- 
ers \={X k ^} and W(b) = Z P (b) is the partition sum of the statistical ensemble associated 
with fixed values b of all the observables. In the following we will call it the extensive 



ensemble using the word "extensive" in the weak sense of "observable" ||14|| , meaning that 
observables are controlled on an event-by-event basis in the extensive ensemble. The two 
partition sums are related through the usual Laplace transform Z x = f db W{b) exp(— A • b). 
Eq. (JTJ) clearly demonstrates that the convexity anomalies of the thermodynamical potential 
logVT can be traced back from log P x {b). Indeed, the equation of states of the extensive 
ensemble can be related to the derivative of log P x (b) 

^,^|M = ^|M +V (2) 

ob k db k 
while the curvature matrix of logp^ is the one of log P x (b) 

Ckk ,(b) = ^S. = dHogWjb) = d 2 log P x (b) ^ 
kk ' db k , db k db k , db k db k , 

It should be noticed that these relations are valid for every set of Lagrange multipliers A 
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because in a finite system the probability distribution covers the entire accessible space of b 
and so contains the whole thermodynamics of the extensive ensemble. The only limitation 
can be a practical one for experiments or simulations: to accumulate enough statistics at 
every location. 

Let us consider the simplified case of a unique direction 6 fc = b and so a unique Lagrange 
multiplier A^ = A. If W (b) has an abnormal curvature, then the EOS A (b) = dblogW (b) 
presents a back-bending. For this extensive ensemble where b is the control parameter, co- 
existence can be defined as the region where A is associated to three values of b because of 
the anomalous curvature. The definition of a first order phase transition as the occurrence 
of a negative heat capacity || |J is readily obtained if one identifies S = log W as the 
microcanonical entropy, b as the energy and A = (3 the inverse canonical temperature. The 
generalization to the EOS of any thermodynamic potential involving an extensive variable b 
H is also obvious. For example, a bimodal grand-canonical number of particle distribution is 
equivalent to a negative chemical susceptibility in the canonical ensemble, a bimodal distri- 
bution of magnetization to a negative magnetic susceptibility in the constant magnetization 
ensemble, a bimodal volume distribution in the isobar ensemble to a negative compressibility 
of the isochore ensemble, ... 

For A in the region of anomalous curvature the associated probability distribution presents 
two maxima and a minimum. In the intensive ensemble eq.(|l|) where A is controlled, coex- 
istence is then signalled by the bimodality of the probability distribution and the value of A 
for which the two maxima have equal height is the first order transition point, analogous to 
the usual Maxwell construction. Therefore, the definitions of a first order phase transition 
through the occurrence of an abnormal curvature of the thermodynamical potential of the 
extensive ensemble, or through the presence of a bimodal event distribution in the associated 
intensive ensemble are strictly equivalent. 

Let us now consider the definition based on the zeroes of Z^, the partition sum of 
the intensive ensemble eq.(^), in the complex A plane. We shall now work out the relations 
between these zeroes and the concavity of P%. Let us continue for simplicity with a single 
observable bk = b and the associated Lagrange multiplier A& = A. Looking at the zeros of Z\ 
in the complex A plane, we first show that a bimodal distribution corresponds to a partition 
sum fulfilling the Yang Lee theorem in the double saddle point approximation |§. Then we 
will demonstrate the reciprocal: a distribution of zeroes which fulfills the Yang Lee theorem 
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is associated to a bimodal probability distribution for b . 

Using equation (fl|) we see that the partition sum for a complex parameter 7 = A + it] is 
nothing but the Laplace transform of the probability distribution P\ (b) for any arbitrary 
parameter Ao JHJ 

Z 7 = J dbZ Xo P Xo (6) e-^- x ^ b = J db PXo (6) e^ b (4) 

where 7 = 7 — Ao and we have defined p x (b) = Z\ P\ {b) . We can also write the partition 
sum at a complex point 7 = A + irj as a Fourier transform of the probability distribution 

Px (&) 

Z y = J db p x (b) e~ irib (5) 

If P\ (b) is monomodal we can use a saddle point approximation around the maximum b\ 
giving Z 7 = e^ 7 ^) , with 

7 (b) = \og Px (b) - i V b + (b) + l - log (2vra 2 (6)) (6) 

where a~ 2 = d%logp\ (b). However, if W\ (b) has a curvature anomaly it exists a range 
of A for which the equation db log W\ = A — Ao has three solutions, one minimum and 
two maxima b^ and b^ . Let us first split the probability distribution into two normal 
components 

Px(b)=p^(b)+p^(b) (7) 
p^ (6) being peaked at b^ . The partition sum reads 

Z 1 = + e 07> = 2e 07 cosh(<^) (8) 

where 20+ = 7 + 7 2 ^ and 20~ = 7 — (jxy . We can now use a double saddle point ap- 



proximation around the two maxima which will be valid close to thermodynamical limit fl5 
and we get </> 7 = 7 (&^) according to eq. (H). The zeros of Z 7 then correspond to 

20; = i(2n + 1)tt (9) 

The imaginary part is given by rj = (2n + l)ir/Ab where Ab = b£ — b^ is the jump in 
b between the two phases. For the real part we should solve the equation dt((j)~) = 0. In 
particular, close to the real axis this equation defines a A which can be taken as Ao- If the 
bimodal structure persists when the number of particles goes to infinity, the loci of zeros 
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corresponds to a line perpendicular to the real axis with a uniform distribution as expected 
for a first order phase transition [|K| . 

Let us now work out the necessary condition and show that a uniform distribution of 
zeroes perpendicular to the real axis with a density linearly increasing with the number of 
particles N implies a bimodal probability distribution. Let us denote zeros as 

7 n = A + *(2n + l)^ (10) 

such that the interval 2ir/8 contains N uniformly distributed zeroes in agreement with the 
unit circle theorem |jI0fl . Since all the zeros of Z 7 are periodically distributed one can define 
an analytic function f(^y) such that 

Z 7 = 2coah(( 7 -Ao)^)/(7) (11) 

Indeed since Z 7 is analytic, f(j) could in principle diverge only on the zeroes of Z 7 ; in these 
points however eq.flTTD shows that / is proportional to <9 7 Z 7 , i.e. it is analytic. On the line 
of zeroes, relation ( |11| ) reduces to Zx +ir] — 2 cos {r]N5/2) f(X Q + ir]). Using the fact that the 
partition sum along the line of zeroes is the Fourier transform of the reduced probability 
distribution at the transition point Ao (eq.(|5|), we can use the inverse Fourier transform to 
get the distribution p Xo (6) at the transition point 

Px (6) = ^ J d V Z Xo+iv (12) 
Equation (|12"D can then be rewritten as 

Pa (6) = 9xo (b + + g Xo (b - ( 13 ) 



9x(b) = ^ / d V f(j) e^ b 



where 

is a distribution. Indeed if we compute p x (b) a little above A , A = A + e we get 

PXq+£ (b) = e^g Xo+£ (& + ^) + e"^/ W (b - ^) (14) 
For large N only the first term survives 

p Xo+£ (b) c e^g Xo+E (& + S {b) (15) 
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which is the first term in the distribution at the transition point eq.(|T2[). If we conversely 
compute p\(b) a little below the transition point we get 

p Ao _ £ (b) * e^g^ £ P ( Ss (P) (16) 

which is the second term in the transition point distribution eq . fll2|) . At the transition point 
px (b) is the sum of two shifted (identical) distributions, p^ and p^\ If g\ is a normal (i.e. 
monomodal) distribution the central limit theorem guarantees that for a large number of 
particles its width will scale as V~N, i.e. will grow slower than the distance between the two 
peaks that scales as N. This implies that (for not too small N) p\ (b) is bimodal. This 
stays true if g\ itself is bimodal (or multimodal), with the only difference that the total 
distribution p\ (b) will then represent the coexistence between more than two phases. Let 
consider the simplest case of a bimodal structure for p Xo (6), close to the thermodynamical 
limit; a little before the transition point only the lowest peak of the b distribution is present, 
while a little above only the second one remains. At the transition both are present and the 
most probable b jumps from b < = b + N5/2 to b > = b — N5/2 where b is the maximum of 
9\ (b) ■ 6 represents the discontinuity of the variable b per particle. If b is the energy 5 is 
the latent heat per particle. 

The important difference between finite systems and the thermodynamic limit is that in 
the latter case the probability distribution is bimodal only at the transition point Ao- On the 
other hand in finite systems, in an interval AA of non-zero measure around the transition 
point A the two phases coexist, i.e. the distribution is bimodal, each peak being associated 
with a phase having a finite probability of occurrence. Eq.flHD shows that this interval is 
the larger the smaller is 5, i.e. the lower is the density of zeroes. This extension of the phase 
transition point to a region in which the intensive ensemble presents a bimodal distribution 
makes a direct measurement of phase coexistence possible [15], contrary to the common 
belief that phase transitions would only be loosely defined out of the thermodynamic limit. 
Moreover, the fact that the distribution is non zero between the two maxima 6< and b > in 



the case of small systems (see eq. fll4|) ) implies that microstates can be accessed that do not 
exist at the thermodynamic limit in the intensive ensemble. These states are specific of the 
coexistence region and can lead to spectacular phenomena as negative heat capacity and 
negative compressibility || in the extensive ensemble. 
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This demonstration can be extended to account for more important finite size effects fl5 
when zeroes 7„ are still periodic but not yet uniformly distributed on a straight line |2|, |T5 . 
Indeed, it is important to remark that the distribution of zeroes is periodic in the imaginary 
direction before the thermodynamic limit. If for example we consider the grancanonical 
ensemble, the Lagrange parameter A = —a = —/3fi represents the logarithm of the fugacity 
z = e a and the observable is the number of particles b = N. The partition sum Z a oc e aN in 
the complex a plane is periodic for any arbitrary value of N. In the more general case, this 
will not be true for very small systems since the extensive variables bk are not in general 
proportional to the number of constituents of the system. However, if the forces are short- 
ranged and iV is sufficiently high, bk oc N, which guarantees the periodicity of Z\ in the 
complex A plane. This periodicity is a constraint on the distribution of zeroes. 

To take into account finite size distortions, we can introduce a transformation 

777.(7) = 1 ~ Am(7) (17) 

which is the identity on the real axis and which maps the zeroes 7 n of the partition sum 
onto a uniform density perpendicular to the real axis 

m(7n) =A + 7(277 + 1)^ (18) 

As before we can introduce an analytic function ^(7) such that 

Z, = 2 cosh ((m (7) - A ) *pj 7( 7 ) (19) 

Using this expression in the inverse Fourier transform eq.([P2"D the probability distribution 
at the transition point reads 

PXo®=&{b+™)+g^(b-™) (20) 

where 

#± (b) =N J drif(X + i v ) e TAm(\ +i V )N6/2 e i v b ^ 

If a phase transition is present at the thermodynamic limit, according to the Yang-Lee theo- 
rem Am goes to zero faster than iV such that the two converge to the same distribution 
gx , and P\ (b) tends to eq. (p~3|) . Following the discussion of eq. (|20|) , this means that a 
deformed distribution of zeroes that converges to a straight line of equally spaced zeroes 
perpendicular to the real axis gives rise to a bimodal probability distribution function. 
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In finite systems, if Am is not zero, the sum (|20|) may or may not present a concavity 
anomaly depending on the actual properties of Am. This is true irrespectively of the fact that 
the straight distribution of zeroes eq.flTOD is asymptotically reached (i.e. a first order phase 
transition exists in the bulk) or not. In this case, the physical phenomenon could be classified 
differently with the zeroes of the partition sum or with the concavity criteria j|, [|. Indeed, 
we have demonstrated their equivalence only for large systems. 

In conclusion, our topological definition of phase transitions through the bimodality of 
the probability distribution function, generalizes the definition based on the microcanonical 
entropy H, || and is equivalent to the definition based on convexity anomalies of generic 
thermodynamic potentials || for any number of particles. We have demonstrated that a 
bimodality is equivalent to the Yang Lee definition of phase transitions close to and at the 
thermodynamical limit. A uniform distribution of zeroes perpendicular to the real axis of a 
Lagrange multiplier A at the position Ao implies that the associated distribution of events 
is bimodal for this Ao , and conversely a bimodal distribution at a given Ao generates a 
uniform distribution of zeroes with a real part Ao and a regularly spaced imaginary part. 
This means that the occurrence of bimodalities can be considered as a valid generalization of 
the concept of phase transition to systems of any size. This topological definition is extremely 
powerful since it gives access to the order parameters defined as the observables for which the 
distribution is bimodal. We have demonstrated that this definition agrees with the definition 
based on the zeroes of Z when the bimodal distribution can be approximated by two 
gaussians, i.e. in the double saddle point approximation. This means that the different 
definitions are coherent for N big enough that fluctuations around the two phases are at the 
gaussian level, but they can differ for smaller systems. The compatibility of the different 
definitions at the thermodynamic limit means that the great effort done by the different 



research groups to recognize and classify phase transitions in finite systems [g, ||, plj is not 
only rigorous but also consistent with macroscopic thermodynamics. 
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